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Abstract 

We investigate the resonant leptogenesis scenario in the minimal B — L extended standard 
model (SM) with the B — L symmetry breaking at the TeV scale. Through detailed analysis of 
the Boltzmann equations, we show how much the resultant baryon asymmetry via leptogenesis 
is enhanced or suppressed, depending on the model parameters, in particular, the neutrino Dirac 
Yukawa couplings and the TeV-scale Majorana masses of heavy degenerate neutrinos. In order 
to consider a realistic case, we impose a simple ansatz for the model parameters and analyze 
the neutrino oscillation parameters and the baryon asymmetry via leptogenesis as a function 
of only a single CP-phase. We find that for a fixed CP-phase all neutrino oscillation data and 
the observed baryon asymmetry of the present universe can be simultaneously reproduced. 
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1 Introduction 



The origin of the baryon asymmetry in the present universe is one of the big mysteries in 
cosmology. The ratio of the baryon (minus anti-baryon) density % to the entropy density s 
has been measured with the precision at 10% level by the WMAP satellite experiment p], 

Y B = — = 0.87 x 10~ 10 . (1) 

s 

It might be the most attractive if the origin of the baryon asymmetry can be explained within 
the context of the SM, the electroweak baryogenesis [2]. In order for this scenario to work, a 
strong first order electroweak phase transition is necessary in the early universe. However, the 
Higgs potential satisfying the current lower bound on the SM Higgs boson mass [3] is not likely 
to show this strong first order phase transition and hence, the SM electroweak baryogenesis is 
almost ruled out. 

An appealing alternative is the leptogenesis scenario which is also intimately related 
with the smallness of the neutrino masses through the seesaw mechanism [5]. A most widely 
accepted scenario is to extend the SM by introducing the right-handed Majorana neutrinos with 
masses around an intermediate scale whose out-of-equilibrium decays create lepton asymmetry 
in the universe. The lepton asymmetry is converted to the baryon asymmetry through the 
(B + L)-violating sphaleron transitions [61 [7] with the conversion rate [8] 

_ 8N f + AN H _ _28 
B ~ 22N f + 13N H L ~ 79 [) 

where we have taken Nf = 3 and Nh = 1 are the numbers of fermion families and Higgs doublets 
in the SM. In normal thermal leptogenesis, there is a lower bound on the mass of Majorana 
neutrinos > 10 10 GeV [9] in order to create sufficient amount of the baryon asymmetry. If it is 
the case, it is hopeless to directly observe the heavy neutrinos at high-energy colliders in the 
near future. 

Many models beyond the SM have been proposed, which may be realized at the TeV scale 
and hence accessible to the Large Hadron Collider (LHC) currently in operation and more fu- 
ture colliders such as the International Linear Collider. Among many models, in this paper, we 
consider the minimal gauged B — L extended SM. This is an elegant and simple extension of the 
SM, in which the right-handed neutrinos of three generations are necessarily introduced for the 
cancellation of the gauge and gravitational anomalies. In addition, the mass of right-handed 
neutrinos arises associated with the U(l) b-l gauge symmetry breaking and the seesaw mech- 
anism is automatically implemented. In the view point of LHC physics, it is very interesting if 
the B — L symmetry breaking scale lies around TeV so that the B — L gauge boson (Z' boson) 
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and the right-handed neutrinos can be discovered in the near future [10]. Recently, we have 
proposed the minimal B — L model with the classical conformal invariance [II] and showed 
that the B — L symmetry breaking in this model is naturally realized at the TeV scale when 
the B — L gauge coupling constant is the same order of magnitude as the size of the SM gauge 
coupling constants [T2] . 

Although the minimal B — L model at TeV is a very attractive scenario, the normal thermal 
leptogenesis scenario cannot work because the mass scale of the right-handed neutrinos is far 
below the bound, 10 10 GeV mentioned above. In this case, the CP-asymmetry parameter, 
which is roughly proportional to Dirac Yukawa coupling squareds is too small to give sufficient 
amount of baryon asymmetry in the universe. However, it has been found that when two 
right-handed neutrinos have almost degenerate masses, there is an enhancement of the CP- 
asymmetry parameter [13], and this enhancement can make the leptogenesis scenario viable 
even if the mass scale of the right-handed neutrinos lie around TeV, the resonant leptogenesis 
[13] . The maximum enhancement is achieved when the mass splitting between two right-handed 
neutrinos is comparable to the decay width of either right-handed neutrinos. By tuning the mass 
splitting between two right-handed neutrinos, even a CP-asymmetry parameter of order unity 
can be obtained in principle. However, it is still non-trivial whether the minimal B — L model 
at the TeV scale can reproduce the observed baryon asymmetry because, as we will discuss 
later in detail, the creation of the lepton asymmetry via decays of right-handed neutrinos is 
highly suppressed in the presence of the \](1)b-l gauge interaction with Z' boson mass at the 
TeV scale [15]. 

In this paper, we investigate in detail the resonant leptogenesis scenario in the minimal B—L 
extended SM with the B — L symmetry breaking at the TeV scale. Through detailed analysis 
of the Boltzmann equations with a variety of model-parameter sets, we show how much the 
resultant baryon asymmetry via leptogenesis is enhanced or suppressed, depending on model 
parameters, in particular, neutrino Dirac Yukawa couplings and TeV-scale Majorana masses of 
heavy degenerate neutrinos. In order to consider a realistic case, we impose a simple ansatz for 
model parameters and analyze the neutrino oscillation parameters and the baryon asymmetry 
via leptogenesis as a function of only a single CP-phase. We find that a fixed CP-phase can 
simultaneously reproduce all neutrino oscillation data and the observed baryon asymmetry in 
the present universe. 

The paper is organized as follows. In the next section, we give a brief review on the minimal 
B — L model and the natural realization of the B — L symmetry breaking at the TeV scale. In 
section 3, we analyze in detail the resonant leptogenesis at the TeV scale by numerically solving 
the Boltzmann equations with various parameter sets. We show how the generated baryon 
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Table 1: Particle content: In addition to the SM particles, right-handed neutrinos N (i = 1, 2, 3 
denotes the generation index) and a complex scalar $ are introduced. 

asymmetry depends on the model parameters such as Dirac Yukawa coupling, right-handed 
neutrino mass spectrum, etc. In section 4, we investigate more realistic parameter choices so 
as to reproduce the neutrino oscillation data. We introduce two right-handed neutrinos and 
a simple ansatz among the parameters, by which the neutrino oscillation parameters and the 
baryon asymmetry are determined by only a single CP-phase. We find that there exists a CP- 
phase which simultaneously reproduces the neutrino oscillation data and the observed baryon 
asymmetry. The last section is devoted for conclusions. Formulas used in our analysis are listed 
in Appendix. 

2 The Minimal B — L Model at TeV 

The minimal B — L extended SM is based on the gauge group SU(3) c xSU(2) L xU(l)y xU(1)b_l 
with the particle contents listed in Table 1. The right-handed neutrinos (iVj) of three generations 
are necessarily introduced by which all the gauge and gravitational anomalies are canceled. 
The SM singlet scalar field ($) works to break the U(1)b-z gauge symmetry by its vacuum 
expectation value (VEV), ($) = Ub-l/v2- Once the B — L gauge symmetry is broken, the Z' 
boson acquires mass, 

% = 2g B _ L v B -L, (3) 

where qb-l is the B — L gauge coupling. The current experimental bound was found to be 

v B -L > 3 TeV US]. 

The Lagrangian relevant for the seesaw mechanism is given by 

C D //X ///I - iyJr^k + h.c, (4) 

where without loss of generality, we work on the basis in which the second term is diagonalized 
and y l N is real and positive. The first term gives the Dirac neutrino mass term after the 
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electroweak symmetry breaking (mp = yi)(H)), while the right-handed neutrino Majorana 
masses are generated through the second term associated with the B — L gauge symmetry 
breaking: 



The B — L symmetry breaking scale is determined by parameters in the Higgs potential and 
in general it can be taken to be any scale as long as the experimental bound Vb-l ^ 3 TeV [T5] 
is satisfied. As discussed in the previous section, we assume the B — L symmetry breaking at 
the TeV scale in this paper, and the masses of Z' boson and right-handed neutrinos lie around 
TeV. In fact, it has been pointed out in [TT| [T2] if we impose the classical conformal symmetry 
on the minimal B — L model, the B — L symmetry breaking can be naturally realized at the TeV 
scale. In the rest of this section, we would like to briefly review the classically conformal B — L 
extended standard model proposed in [11]. However, since the classical conformal invariance 
is not important for the leptogenesis scenario (except that it naturally leads to the TeV scale), 
readers can skip to the next section. 

We first note that because of its chiral nature, the SM Lagrangian at the classical level 
possesses the conformal invariance except for the Higgs mass term, closely related to the gauge 
hierarchy problem. Bardeen has argued [TT] that once the classical conformal invariance and 
its minimal violation by quantum anomalies are imposed on the SM, it could be free from the 
quadratic divergences and thus the gauge hierarchy problem. If the mechanism really works, 
we can directly interpolate the electroweak scale and the Planck scale. Since the classical 
conformal symmetry forbids the mass term in the Higgs potential, the electroweak symmetry 
should be broken radiatively through the Coleman- Weinberg (CW) mechanism [18]. Although 
this is an attractive scenario, the effective Higgs potential is found to be unbounded from below 
because of the large top Yukawa coupling and therefore the classically conformal SM cannot be 
a realistic scenario. 

In [TT], we proposed a classically conformal minimal B — L model and showed that the 
B — L gauge symmetry breaking is successfully achieved via the CW mechanism and then, 
this breaking triggers the electroweak symmetry breaking. Because of the CW mechanism, the 
SM singlet Higgs boson associated with the B — L symmetry breaking is much lighter than Z' 
boson, 



where cxb-l = 9b-lI '(^tt), and a % N = (y l N ) 2 / (4n) . This formula also indicates the upper bound 
on a l N to keep the vacuum stability, m 2 , > 0. Assuming a hierarchical Majorana mass spectrum, 




(5) 
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for example, we find the upper bound on the heaviest right-handed neutrino mass as 



E4 3 4 



(7) 



Once the B — L symmetry is broken, the Z' boson and the right-handed neutrinos acquire 
masses at the B — L symmetry breaking scale. Their masses contribute to the effective mass of 
the SM Higgs doublet through quantum corrections. The naturalness argument, namely, the 
quantum corrections should not exceed the electroweak scale so far, leads to an upper bound 
on the B — L symmetry breaking scale. Two-loop corrections with Z' boson and top quarks 
are found to be dominant, and we have concluded [12] that the B — L symmetry breaking scale 
should be around TeV. 

3 Leptogenesis in the minimal B — L model at TeV 

Now we study baryogenesis via leptogenesis in the minimal B — L model at TeV. The lep- 
ton asymmetry in the universe is generated by the CP-violating out-of-equilibrium decays of 
right-handed neutrinos, and this asymmetry is converted to the baryon asymmetry through the 
sphaleron process with the conversion rate Y# = — (79/28)Y^. The generated baryon asym- 
metry is evaluated by solving the Boltzmann equations. When the Majorana masses of three 
right-handed neutrinos are largely different as usually assumed, it is sufficient to consider the 
Boltzmann equations only for the lightest right-handed neutrino. This is because the lepton 
asymmetry generated by heavier right-handed neutrinos are washed out by the inverse-decay of 
the lightest right-handed neutrinos before its out-of-equilibrium decay [19]. However, in the res- 
onant leptogenesis scenario, (at least) two right-handed neutrinos are degenerate in mass and it 
is generally not clear whether analysis of the Boltzmann equations with only one right-handed 
neutrino is sufficient. As we will see later, it can be essential for general cases to consider the 
Boltzmann equations for multiple right-handed neutrinos. 

We begin our analysis with the Boltzmann equations in one-flavor approximation^, 
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Throughout the paper, our notation follows Ref. Q]5] 
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Figure 1: Feynman diagrams corresponding to each 7s. 
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Figure 2: Right-handed neutrino decay at tree and one-loop levels. 



where is the yield (the ratio of the number density to the entropy density s) of the (lightest) 
right-handed neutrino, Y^ q is the yield in thermal equilibrium, temperature of the universe 
is normalized by the mass of the right-handed neutrino z = Mi/T, H(Mi) is the Hubble 
parameter at T = Mi, ei is the CP-asymmetry parameter, and 7s are the space-time densities 
of the scatterings in thermal equilibrium. Diagrams in Fig. [1] show processes corresponding 
to different 7s in the Boltzmann equations, whose explicit forms are listed in Appendix. As 
a good approximation, we neglect masses for all particles involved in the processes, except for 
the right-handed neutrinos and the Z' boson having the TeV-scale masses. The yield of the 
right-handed neutrino obeys the first equation, while the send equation determines the B — L 
number created by the out-of-equilibrium decay of the right-handed neutrinos with a non-zero 
CP-asymmetry parameter. In our numerical studies with input parameters given below, we 
can check that only 7^ and jz> among the space-time densities have important effects on the 
final results while the others are negligible. 

The CP-asymmetry parameter associated with the decay of right-handed neutrino JVj is 
defined as 



£\ [r -> £jH) -r(N^ efH*)] 



(9) 



£ . [r £jH) + r (Ni -> efH*)] ' 

which is generated by the interference between the tree and one-loop diagrams shown in Fig. [2j 
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The general formula is given by [131 E 
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corresponds to the vertex correction of the second diagram in Fig. [2] while 
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11] 



(12) 



(am 2 ) +M 2 r 2 

is from the self-energy corrections of the third diagram in Fig. [2] with the decay width and the 



mass difference defined as 
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(13) 



In the terminology of K° — K° mixing, the first contribution is the so-called direct CP violation 
while the second one the indirect CP violation. The indirect CP violation occurs because the 
mass eigenstates and the CP eigenstates are generally different. 

When the right-handed neutrinos have a hierarchical mass spectrum (Mi -C M 2j 3), the con- 
tributions from Vj and Sj are comparable and the CP-asymmetry parameter is approximately 
given by [20j 



f 7 ! 



3 m v M\ 



sin S ~ 10" 



m,. 



Mi 



0.05eV/ V10 10 GeV 



sin 5, 



(14) 



167T V' 

where m v is the light neutrino mass eigenvalue, and 5 is the CP-phase and assumed to be of 
order one. The baryon asymmetry of the universe is parameterized as 



Y B = k , 
9* 



(15) 



where g* = (9(100) is the number of relativistic degrees of freedom in the early universe, and 
k is the efficiency factor determined by solving the Boltzmann equations, independently of the 
CP-asymmetry parameter. In the leptogenesis scenario without the Z' gauge boson, this factor 
is roughly estimated as [9] 



k ~ 2 x 10" 



0.05eV 



1.1 



m,. 



(16) 



Using these estimations, we arrive at the conclusion that the Majorana mass Mi has to be larger 
than 10 10 GeV to give the observed value Yb ~ 10 -10 . However, note that this conclusion 



7 



is based on the formula of the CP-asymmetry parameter in the case with the hierarchical 
right-handed neutrino mass spectrum. In fact, when two right-handed neutrinos are almost 
degenerate, the CP-asymmetry parameter can be enhanced. 

Now, suppose that two right-handed neutrinos are almost degenerate. In this case, it is 
easy to note that there is a parameter region which can dramatically enhance Sj. Maximum 
enhancement occurs for AM^ ~ MiTj <C Mf, so that Sj ~ MjjVj 3> 1. In this case, the 
CP-asymmetry parameter is given by 



Im 










33 



(17) 



which can, in principle, be of order unity. The leptogenesis scenario with this enhancement 
of the CP-asymmetry parameter is called the resonant leptogenesis [TJ]. This enhancement 
is crucial to realize the observed baryon asymmetry when the right-handed neutrino mass is 
significantly smaller than 10 10 GeV, such as the TeV scale which is of our main concern in this 
paper. 

Analysis of the Boltzmann equations in one-flavor approximation 

Now we analyze the Boltzmann equations in Eq. (jSJ) to see how much baryon asymmetry can 
be generated in the B — L model at the TeV scale. In order to understand the response between 
the model-parameters involved in this analysis and the resultant baryon asymmetry, we first 
consider a model in one-flavor approximation with the parameterization of the decay width as 

Ti = V ^M X (18) 

with a real free parameter yjj, while the other parameters are fixed as follows: 

Cl = 0.01, ol b~l = 0.006, m z > = 3 TeV, M t = 2 TeV. (19) 

Then, we numerically solve the Boltzmann equations with the boundary conditions 

Y Nl (0) = Y%{0), F B _ L (0)=0. (20) 

The lepton asymmetry generated by the right-handed neutrino decays is converted into the 
baryon asymmetry via the sphaleron process while the process is in thermal equilibrium. In 
our analysis throughout the paper, we evaluate the resultant baryon number at the freeze-out 
temperature of the sphaleron process, T sp h — 150 GeV [21], where the conversion of the lepton 
number to the baryon number is terminated: 

Y B = ^Y B _ L (z sph ), (21) 
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log yD 2 

Figure 3: The baryon asymmetry generated in the minimal B — L model (solid red line) as 
a function of the Dirac Yukawa coupling. The dotted (blue) line corresponds to the result of 
baryon asymmetry in the absence of the B — L gauge interaction, while the dashed (green) line 
is for the case with jN,t,& — as well as as-L = 0. 

where Y B _ L (z sph ) is the numerical solution of the Boltzmann equations at z sp h = Mi/T sph . 

The resultant baryon asymmetry is depicted in Fig. 3 as a function of yp (solid line). For 
comparison, results for the cases with as-i — (dotted blue line) and with ^N,t,<5> — as 
well as as-L — (dashed green line) are also shown. For a small y 2 D < 10~ ia5 , we can see 
that the generation of the baryon asymmetry is suppressed in the presence of the Z' boson and 
lN,t,<& processes. Although the suppression by the Z' boson effect dominates, the ■jN,t,<s> process 
mediated by the Majorana Yukawa coupling (y N ) causes a dramatic reduction in generating the 
baryon asymmetry even in the absence of the Z' boson effect. In the region, Yb is growing as 
yoi and a larger yo generates a larger baryon asymmetry against the Z' boson and Majorana 
Yukawa coupling effects. On the other hand, for y 2 D > 10~ 10 , the effect by Dirac Yukawa 
coupling dominates over the Z' boson and Majorana Yukawa coupling effects, and all lines 
become well-overlapping. In this region, however, Yb is suppressed by the washing-out process 
via the inverse-decay process. In the dashed (green) line, Yb becomes smaller as y 2 D is lowered 
f° r Vd ~ 10~ 15 GeV, nevertheless •Jnj,® = and as-L = 0. This is because the generation of 
lepton number is too slow with such a small Dirac Yukawa coupling, and the sphaleron process 
freezes out before the completion of the whole lepton number generation. 

Fig. H] shows the results for different values of M\ as a function of y 2 D while the other 
parameters are kept the same. We can see, for M\ > 1.5 TeV, a similar behavior to the 
result shown in Fig. 3. For a relatively small Mi < 1.5 TeV, the resultant Yb becomes almost 
independent of y 2 D even for a larger Yp. This is because the freeze-out of the sphaleron process 
occurs and thus the conversion of the lepton number to the baryon number is terminated before 
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Figure 4: The baryon asymmetry as a function of Yukawa coupling for different values of the 
right-handed neutrino mass, Mi = 1 TeV (solid red line), Mi = 1.5TeV (dashed green line) , 
Mi = 2 TeV (dotted blue line) and Mi = 2.5 TeV (solid black line). 

the washing-out process becomes effective. 
2-flavor analysis 

We have so far investigated the leptogenesis by solving the Boltzmann equations with one- 
flavor right-handed neutrino. This treatment is justified when the right-handed neutrinos have 
a hierarchical mass spectrum. In the resonant leptogenesis, two right-handed neutrinos are 
degenerated in mass and thus, it is non-trivial whether one-flavor analysis is actually a good 
approximation or not. Here, we generalize our analysis to the 2-flavor case and clarify when 
the 1-flavor analysis is justified. 

As we have discussed in the previous section, the CP-asymmetry parameter e« is maximally 
enhanced when two right-handed neutrinos are almost degenerate and their mass squared dif- 
ference is AMfj ~ MiTj. Thus, in our analysis for the two-flavor case, we set the mass difference 
as AMf 2 = M 1 T 2 . Note that since M li2 > I\ 2 and thus Mi ~ M 2 , the case with AMf 2 = M^i 
is essentially the same as the case with the exchange 1 -H- 2. The CP- asymmetry parameters 
are given by 
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(22) 



where we have used the relations (hdUd)^ = (VdUd)^ an< ^ AM^ 2 = — AM| 3 
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Figure 5: The ratio Fjl^ avor /F^I^ avor as the Dirac Yukawa coupling squared y 2 Dl for different 
choices of ^ = 10 (red solid line) j£ = 100 (green dashed line) and j£ = 1000 (blue dotted 
line) . 



We consider the following three cases; 

1. r\ >r 2 

The CP asymmetry parameter e 2 — 2eiP 2 /Ti <C e 1; and the baryon asymmetry is gen- 
erated dominantly by the N\ decay. In addition, the washing-out process by the inverse 
decay of iV 2 is also negligible to that by N\. Therefore, analysis with only one- flavor 
right-handed neutrino JVi is sufficient in evaluating the resultant baryon asymmetry. 

2. r\ ~ r 2 

Clearly, two right-handed neutrinos are almost identical, so that one-flavor analysis is 
sufficient, but the resultant baryon asymmetry should be twice of that obtained in one- 
flavor case. 

3. r x < r 2 

The CP asymmetry parameter e 2 ~ 2ei^ <C e\ and hence the generation of baryon 
asymmetry by iV 2 decays are negligible. However, the washing-out effect by the inverse 
decay of iV 2 can be efficient and the generated baryon asymmetry can be drastically 
reduced, depending on the value of T 2 . 

Therefore, for the cases 1 and 2, one-flavor analysis is sufficient to evaluate generated baryon 
asymmetries, while the case 3 is non-trivial and we need to solve the Boltzmann equations with 
two-flavor right-handed neutrinos. 
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The Boltzmann equations with two flavor right-handed neutrinos are given by 

(l-0- + ((iH^ 

Here we have omitted 7s mediated by the Higgs boson and the right-handed neutrinos, because 
these effects are, in fact, negligible. However, all processes in Fig. 1 are taken into account in our 
numerical analysis. With the initial conditions, Y^O) = Y^(0) (i = 1, 2) and Yb-l(Q) = 0, we 
solve these equations and evaluate the baryon number at T sp h = 150 GeV. In order to compare 
results to those in the one-flavor case, we parameterize the decay width as Ti = y^Mi/ 
corresponding to the one-flavor case. For fixed values of T 2 /Ti, we show the ratio of the baryon 
asymmetry to the one obtained in the one- flavor analysis, Fjlf vor /l^:f vor , in Fig. El For 
a small yp, we can see that the one-flavor analysis is a good approximation. On the other 
hand, the washing-out process by the inverse-decay of N 2 is very effective for a large yo, and 
the baryon asymmetry is very much suppressed than the result obtained in the one-flavor 
analysis. As is expected, the baryon asymmetry is more suppressed as the ratio T2/T1 becomes 
larger. Therefore, in order to obtain the correct result for baryon asymmetry via the resonant 
leptogenesis, analysis with two (or more) flavors can be essential in the general case, especially, 
when 6i <C 6j but Tj 3> Tj for two almost degenerate right-handed neutrinos, Ni and Nj. 
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4 Neutrino oscillation data and resonant leptogenesis 

In the previous section, we have analyzed the resonant leptogenesis in the minimal B — L model 
and investigated the response of the resultant baryon asymmetry to model-parameters such as 
the Dirac Yukawa couplings and the right-handed neutrino masses. It is clearly more interesting 
to consider a realistic model (in other words, a realistic parameterization) which can account 
for the observed neutrino oscillation phenomena. For this purpose, we consider a strategy first 
proposed in [23] in this section. In our analysis, we adopt the current neutrino oscillation data 
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in 2-er range 

7.25 x 10~ 5 < Am 2 2 (eV 2 ) < 8.11 x 1(T 5 , 



2.18 x 10" 3 < |Am 2 3 |(eV 2 ) < 2.64 x 1(T 3 , 
0.27 < sin 2 # 12 < 0.35, 
0.39 < sin 2 #23 < 0.63, 

sin 2 9 13 < 0.040, (24) 

for the standard parameterization of the mixing matrix, 

cos 6>i2 cos 0i3 sin 0i2 cos 0i3 sin0i3e — ' s 

UpMNS = I — sin 0i2— cos 0i2 sin 023 cos 0i3e l5 cos0i2 cos023— sin0i2 sin023 sin0i3e l * sin023COS0i3 

sin 0i2 sin 023— cos 0i2 cos 023 sin 9i^e' s — cos 0i2 cos 023— sin 0i2 cos 023 sin 0i3e 8<5 cos 023 cos 013 

xdiag(e^,e^,l), (25) 

with the Dirac phase 5 and the Majorana phases a\. 

We consider the so-called minimal seesaw [25] in the context of the minimal B—L model, and 
assume that only two right-handed neutrinos are relevant for the neutrino oscillation phenomena 
and leptogenesis. The third right-handed neutrino is assumed to decouple from the neutrino 
oscillation phenomena by some reason. A simple idea is to introduce a discrete Z 2 symmetry 
under which the third right-handed neutrino is assigned to be odd while all the other particles 
in the B — L model even. In the context of the minimal B — L model with the Z 2 symmetry, 
it has been shown [26] that the third right-handed neutrino can be a suitable candidate for the 
cold dark matter with the relic density consistent with observations. 

In the minimal seesaw model, we parameterize the 2x3 Dirac neutrino mass matrix, without 
loss of generality, as 

m D = [ , (26) 

V °4 as «6 / 

where a, and <ftj are real parameters, and we have worked in the basis where both the charged 
lepton mass matrix and the right-handed neutrino mass matrix are diagonalized with real and 
positive eigenvalues. We parameterize the Majorana mass matrix of the two right-handed 
neutrinos as 

M »=( M M l( lV))' (2?) 

where the parameter r should be very small, for example, r ~ T\/Mi or T2/M1 in order to 
realize the enhancement of the CP-asymmetry parameter. Although r is crucial for the resonant 
leptogenesis, such a small r is negligible in fitting for the neutrino oscillation data. 
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Figure 6: The neutrino oscillation parameters, Amu (left panel) and Ami3 (right panel), as a 
function of the CP-phase 03. The observed data in 2-er range are indicated by two horizontal 
lines. 
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The neutrino oscillation parameters, sin 2 # 12 (left panel) and sin 2 # 2 3 (right panel), as 
of the CP-phase 03. The observed data in 2-a range are indicated by two horizontal 



For simplicity, we fix <pi = 2 = in our analysis and introduce an ansatz 
neutrino mass matrix after the seesaw mechanism [5], 

1V1\ 

is diagonalized by the so-called tri-bimaximal mixing matrix 

/ ./I ,/i \ 

V 3 V 3 

Utb = 



\ 



2 / 



that the light 
(28) 



(29) 



in the CP invariant case (03 = 0). As is well-known, this tri-bimaximal mixing matrix gives 
almost the best fit in the oscillation data. 
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Figure 9: The CP asymmetry parameters, e\ (right) and e 2 (left), as a function of r. 

Note that in the minimal seesaw, the rank of the light neutrino mass matrix is two and the 
lightest mass eigenvalue is 0. In the following we consider two cases for the light neutrino mass 
spectrum, namely, the normal hierarchical (NH) case and inverted hierarchical (IH) case. Let 
us first consider the NH case, where we have 

AT = diag(0,m 2 ra ,mH (30) 

with mg H = y/A 

?7i^2 and — | A?n^ 3 1 . According to our ansatz, we first find a solution 
to m v = UtbD^ h U^ b in the CP-invariant case. Among several solutions, we choose, as an 
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Figure 10: The resultant baryon asymmetry in the universe as a function of the CP-phase 03. 
The observed value Yg = 0.87 x 10~ 10 is depicted as the horizontal line. 



For the input values, we use m^ 11 = 8.75 x 10 -3 eV and m^ H = 4.90 x 10~ 2 eV. In this section, 
we fix other parameters as (Xb-l — 0.006, mz< = 3 TeV and M\ = 2 TeV. 

Now we turn the CP-phase 03 on. With fixed a« and M±, the light neutrino mass matrix 
is given as a function of the single parameter 3 [23] (in the approximation with r = 0). In 
Figs. 6—8, the neutrino oscillation parameters are depicted as a function of the CP-phase 03. 
As we expect, the outputs of the oscillation parameters deviate from the values at 03 = as 
the CP-phase is changed, and eventually some of outputs are found to be outside of 2-a range 
of the experimental data. We find the bound on the CP-phase as |03| < 0.5. 

For the resonant leptogenesis, both 3 ^ and r ^ are crucial. Fig. 9 shows the CP- 
asymmetry parameters (ei and 62) as a function of r with 3 = 0.5 for example. For the right 
(left) curve corresponding to t\ (e 2 ), a peak appears around r = 10~ 13 (r = 10~ 14 ). Choosing 
r = 10~ 14 for example, the CP-asymmetry parameters are also given as a function of only 
03. Therefore, we have correlations between neutrino oscillation parameters and the baryon 
asymmetry generated by the resonant leptogenesis through the CP-phase 03 [23] . Interestingly, 
the amount of the generated baryon asymmetry becomes larger as 03 goes away from zero, 
while a large displacement of 03 from zero results in the oscillation parameters inconsistent 
with the experimental data. 



example, 





0. 



(31) 
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Numerical solution of the Boltzmann equations with two flavors in Eq. (1231) is shown Fig. 10 
as a function of 3 . We find that the observed baryon asymmetry Yg = 0.87 x 10~ 10 in the 
present universe is obtained for 3 = 0.35, for which the neutrino oscillation parameters are 
fixed as 



Am 2 2 (eV 2 ) 


= 7.39 x 10 


Am 2 3 (eV 2 ) 


= 2.39 x 10 


sin 2 6*12 


= 0.34, 


sin 2 #23 


= 0.51, 


sin 2 0i3 


= 0.00016. 



-3 



(32) 

They are all consistent with observations. Although a non- vanishing sin 2 0i3 is predicted, it is 
quite small, far below the current upper bound. 

Next we consider the IH case, where the light neutrino mass matrix is diagonalized as 

Z^ = diag(m^,m^,0), (33) 

with m[ H = a/| Am 2 3 | and = \/ Am^ + |Amf 3 |. In the CP invariant case, we choose a 
solution to m v = UtbDI h U^ b as 



d\ = CL2 = 3 ~~ 



2M 1 m{ 



IH 



M m 1JH 

a 5 = a 6 = -W . (34) 

For the input values of m[ H = 4.90 x 10 -2 eV and m 2 = 4.98 x 10~ 2 eV, the neutrino oscillation 
parameters are depicted in Figs. [TTHTBI as a function of 3 . A CP-phase |0 3 | < 0.1 results in 
the outputs of the neutrino oscillation parameters consistent with the experimental data in 2-a 
range. 

Fig. [T4] shows the CP-asymmetry parameters, e\ and function of r for 3 = 0.1. 

Two curves are well-overlapped and the peak appears around r ~ 10~ 13 . Then, we solve the 
Boltzmann equations for r = 10~ 13 , and show the results in Fig. [15j Although the observed 
baryon asymmetry in the universe is generated for 3 = 0.43 > 0.1, the neutrino oscillation 
parameters corresponding to the CP-phase are outside of 2-a range. Therefore, in the present 
scheme, the IH case cannot reproduce the neutrino oscillation data and the observed baryon 
asymmetry simultaneously. 
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Figure 11: The neutrino oscillation parameters, Ami2 (left panel) and Ami3 (right panel), as 
a function of CP-phase 03. The observed data in 2-cr range are indicated by two horizontal 
lines. 




-3-2-10123 -3-2-10123 





Figure 12: The neutrino oscillation parameters, sin 2 9 12 (left panel) and sin 2 #23 (right panel), 
as a function of CP-phase 03. The observed data in 2-cr range are indicated by two horizontal 
lines. 

5 Conclusions 

In this paper, we have investigated a possibility to explain the baryon asymmetry of the universe 
as well as the neutrino oscillation data in the TeV scale B — L model. In the model, the lepton 
asymmetry is generated in the early universe via out-of-equilibrium decays of the right-handed 
neutrinos with the CP-asymmetry parameter, and converted into the baryon asymmetry via the 
sphaleron process. When the mass scale of the right-handed neutrinos is low < 10 10 GeV, the 
enhancement of the CP-asymmetry parameter is crucial in order to generate sufficient amount 
of baryon asymmetry in the universe. The enhancement is realized when two right-handed 
neutrinos are almost degenerated and in this case, the CP-asymmetry parameter can be in 
principle order unity. This scenario is called the resonant leptogenesis. However, it is still 
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Figure 13: sin 2 #13 as a function of the CP-phase 03. 
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Figure 14: The CP asymmetry parameters, e\ and as a function of r. Two curves are 
well-overlapped. 

non-trivial whether the resonant leptogenesis can realize the observed baryon asymmetry in 
the context of the minimal B — L model, because the B — L interaction mediated by the Z' 
boson can dramatically reduce the generation of baryon asymmetry. 

We numerically solved the Boltzmann equations for the resonant leptogenesis in the minimal 
B — L model, and figured out the response between the generated baryon asymmetry and the 
model-parameters such as the neutrino Dirac Yukawa couplings and the right-handed neutrino 
masses. We first analyzed the Boltzmann equations with only one-flavor right-handed neutrino 
and a fixed CP-asymmetry parameter. When the neutrino Dirac Yukawa coupling is small, 
Ud % 10~ ia5 , the amount of the baryon asymmetry becomes larger as the Dirac Yukawa 
coupling is raised. In this parameter region, the Z' boson and Majorana Yukawa coupling 
effects dramatically suppress the generation of baryon asymmetry. For a large Dirac Yukawa 
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Figure 15: The resultant baryon asymmetry as a function of the CP-phase 03, along with the 
observed value Yg = 0.87 x 10~ 10 (horizontal line). 

coupling, the baryon number generation by the right-handed neutrino decay dominates over 
the suppression by the Z 1 boson effect. However, a too large Dirac Yukawa coupling in turn 
suppresses the generation of the baryon asymmetry by the washing-out effect via the inverse- 
decay process. Next, we have analyzed the Boltzmann equations with two-flavor right-handed 
neutrinos and shown that two-flavor analysis can be essential in general cases. With these 
analyses, we have shown that in some areas of the parameter space a sufficient amount of the 
baryon asymmetry can be generated though the resonant leptogenesis in the TeV-scale B — L 
model. 

Finally, we have checked whether these parameters are consistent with the current neutrino 
oscillation data. We have introduced a simple ansatz for the neutrino mass matrices, by which 
the neutrino Dirac Yukawa couplings are determined as a function of a single CP-phase. For 
both the normal hierarchical and inverted-hierarchical mass spectra of the light neutrinos, we 
have shown the correlations between the neutrino oscillation parameters and the generated 
baryon asymmetry via the resonant leptogenesis. In our analysis with the ansatz, a fixed 
CP-phase can reproduce simultaneously the neutrino oscillation data and the observed baryon 
asymmetry in the normal hierarchical case. On the other hand, we cannot find such a CP-phase 
in the inverted hierarchical case. 
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Appendix 

The number density of a particle ip (a right-handed neutrino in our case) with mass 
in the early universe is evaluated by solving the Boltzmann equation of the form [22] . 



dz 



Y^Ya 



v eq v eq 

z 4> Ia ' 



-7 e9 (^ + a + >i+ ] + ■■■) 



YY 

l eq {i+3 + >ip + a + 



yeqyeq 
i 3 



(35) 



where Y^ = n^/s is the ratio of and the entropy density s, z — ^r, and H(m^) is the Hubble 
parameter at a temperature T = m^. The right hand side of Eq. ( 135]) describes the interactions 
that change number of ip, and 7 e,? is the space-time density of scatterings in thermal equilibrium. 
For a dilute gas we take into account decays, two-particle scatterings and the corresponding 
back reactions. One finds, for a decay the particle ifi, 

r*l(*)p 



7D = 7 e9 -M+j + •••) 



n 



eg - 



(36) 



where K\ and are the modified Bessel functions, and r rs is the decay width. For two body 
scattering one has 

T 



7 (ip + a -H- i + j + ■ ■ 



( \fs 

dsa{s)y/sK 1 — 
64vr 4 J {m ^ +ma) 2 V T 



(37) 



where s is the squared center-of-mass energy and the reduced cross section a for the process 
ip + a-^i + j + -- - is related to the usual total cross section a(s) by 



(s 



<T(S 



(3* 



In the following we list the explicit forms of the reduced cross sections used in our analysis 
19] . The reduced cross section corresponding to 7^ is given by 



or 



2vr 1 



a N {s 



sin 4 OM^x 



11 



2x 



x 



X 



+ 



X 



DAx) 



D 1 (x) 2D\{x) 



In (a; + 1) 



where x 



and 



Di(x) — x — 1 
while the one corresponding to jNt is 



with c 



x 



TTliv 



o"v,H s 



mr)m D 



u \2(x+l) 1 x+2 



+ -^ln(x + l' 



(39) 



(40) 



(41) 
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The reduced cross section for the t-channel (and u-channel) process N + iV — > $ + $ mediated 
by the right-handed neutrino is given by 

y%x — 4 / x x 2 — 8x + 16 , x — J xix — 4) \ , , 

8ir x y 2 xy/x(x-4:) x + ^/x(x-A)J 

The reduced cross section for the s-channel process N + l — > i+q mediated by the Higgs doublet 
is given by 

A , . 3i:a 2 m 2 / + \ /x-l\ 2 , . 

^'^MPTiH),!— J • (43) 

while for the t-channel process 



A . . 3ira 2 m 2 / t \ |~x — 1 1 / x — 1 + ?/\ 
a h,t\ s ) = 774 .4. ( m D m^ ) + -In 



(44) 



with y' = The total reduced cross section for the process / + /,$ + $ — T-iV + iV mediated 
by the Z' boson ( / denotes the SM fermions) is given by 



A , 1047T 2 v/iC 3 



where y = and c = flVM^il with the decay width of Z' boson 

m N \ J 

(46) 



~ a B _ L m z , 
1 Z' — 7. 



3( 1--) 0(2,-4) + 13 



y 
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